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Abstract 

We  consider  systems  governed  by  partial  differential  equations  with  spatially  periodic  coefficients  over  un¬ 
bounded  domains.  These  spatially  periodic  systems  are  considered  as  perturbations  of  spatially  invariant  ones,  and 
we  develop  perturbation  methods  to  study  their  stability  and  system  norm.  The  operator  Lyapunov  equations 
characterizing  the  norm  are  studied  using  a  special  frequency  representation,  and  formulae  are  given  for  the 
perturbation  expansion  of  their  solution.  The  structure  of  these  equations  allows  for  a  recursive  method  for  solving 
for  the  expansion  terms.  Our  analysis  provides  conditions  that  capture  possible  resonances  between  the  periodic 
coefficients  and  the  spatially  invariant  part  of  the  system.  These  conditions  can  be  regarded  as  useful  guidelines 
when  spatially  periodic  coefficients  are  to  be  designed  to  increase/decrease  the  norm  of  a  spatially  distributed 
system.  The  developed  perturbation  framework  also  gives  simple  conditions  for  checking  exponential  stability. 


I.  Introduction 

The  terms  distributed-parameter  and  infinite-dimensional  are  used  to  describe  those  systems  in  which  the  state 
belongs  to  an  infinite  dimensional  vector  space  [1].  Such  systems  include,  but  are  not  limited  to,  time-delay 
(retarded)  and  spatially  distributed  systems  [2],  The  latter  includes  systems  in  which  the  dynamics  are  governed 
by  Partial  Differential  Equations  (PDEs)  and  it  is  a  subclass  of  these  systems  that  will  be  the  subject  of  this  study. 
More  specifically,  we  will  analyze  certain  properties  of  spatially  distributed  systems  in  which  the  underlying  PDEs 
have  spatially  periodic  coefficients.  We  refer  to  such  systems  as  spatially  periodic.  Spatially  periodic  systems  have 
many  real  life  applications,  for  example,  in  boundary  layer  and  channel  flow  problems  with  corrugated  walls  and 
in  nonlinear  optics. 

Our  motivation  for  this  work  is  to  study  the  effect  of  spatially  periodic  coefficients  on  stability  and  system 
norms  of  spatially  distributed  systems.  This  can  be  thought  of  as  using  the  periodic  coefficients  as  static  feedback 
controls  for  spatially  distributed  systems.  Eor  example,  in  flow  control  problems  where  one  introduces  corrugated 
wall  geometries  or  spatially  periodic  body  forces,  the  PDEs  that  describe  the  resulting  flow  dynamics  have  periodic 
coefficients  that  are  related  to  either  the  wall  shapes  or  the  spatially  distributed  body  forces.  An  important  objective 
is  to  “design”  such  wall  shapes  or  body  forces  to  obtain  certain  stability  or  instability  properties  of  the  resulting 
dynamics.  There  are  currently  no  systematic  methods  for  achieving  this. 

An  analogy  can  be  made  between  the  present  work  and  the  use  of  time-periodic  coefficients  in  Ordinary 
Differential  Equations  (ODEs).  It  is  known  that  the  introduction  of  time-periodic  coefficients  in  ODEs  with 
constant  coefficients  can  change  the  stability  properties  of  the  Linear  Time  Invariant  (LTI)  system  described  by 
the  original  ODE.  A  useful  picture  is  to  think  of  an  ODE  with  periodic  coefficients  as  an  LTI  system  modified 
by  time-periodic  (memoryless)  feedback.  It  is  known  that  certain  unstable  LTI  systems  can  be  stabilized  by  being 
put  in  feedback  with  periodic  gains  of  properly  designed  amplitudes  and  frequencies  [3].  This  can  be  roughly 
considered  as  an  example  of  “vibrational  control”  [4].  On  the  other  hand,  certain  stable  or  neutrally  stable  LTI 
systems  can  be  destabilized  by  periodic  feedback  gains.  This  is  sometimes  referred  to  as  “parametric  resonance” 
in  the  dynamical  systems  literature  [3].  In  the  above  scenarios,  the  stabilization/destabilization  process  depends  in 
subtle  ways  on  “resonances”  between  the  natural  modes  of  the  LTI  subsystem  and  the  frequency  and  amplitude 
of  the  periodic  feedback.  Although  Eloquet  analysis  can  be  used  to  ascertain  stability  of  the  resulting  periodic 
systems,  it  is  cumbersome  to  use  for  designing  the  periodic  coefficients.  This  requires  an  exhaustive  search 
over  frequencies  and  amplitudes  of  the  periodic  coefficients.  Alternatively,  simple  resonance  conditions  can  be 
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derived  using  a  perturbation  analysis  [3],  which  in  turn  can  be  used  for  a  preliminary  selection  of  the  coefficient’s 
frequency.  In  this  way,  perturbation  analysis  serves  as  a  useful  design  tool. 

In  related  recent  work  [5]  we  developed  computational  tools  to  study  stability  and  system  norms  for  spatially 
periodic  systems.  However,  for  problems  where  the  spatially  periodic  coefficients  are  to  be  designed,  using  these 
tools  involves  a  computationally  expensive  search  over  spatial  frequencies  and  amplitudes  of  the  coefficients. 
Therefore,  our  aim  in  the  present  work  is  to  develop  a  perturbation  analysis  that  can  be  used  to  derive  resonance 
conditions  and  provide  a  useful  design  tool  in  a  similar  manner  to  the  case  of  ODEs  discussed  earlier.  These 
resonance  conditions  can  then  identify  candidate  spatial  frequencies  to  be  used  for  the  periodic  coefficients.  The 
exact  behavior  with  respect  to  amplitudes  can  then  be  ascertained  using  the  full  analysis  of  [5].  In  this  manner 
we  reduce  the  dimension  of  the  search  space  required  for  design  problems. 

Another  challenging  problem  is  checking  the  stability  of  a  spatially  periodic,  or  in  general,  any  infinite¬ 
dimensional  system.  It  is  well-known  that  for  a  finite-dimensional  LTI  system,  the  spectrum  of  the  infinitesimal 
generator  (i.e.,  the  A-matrix)  being  contained  in  the  open  left  half  of  the  complex  plane  is  equivalent  to  exponential 
stability.  In  this  sense  the  spectrum  of  the  infinitesimal  generator  determines  stability.  Therefore  it  is  said  that 
the  system  satisfies  the  Spectrum  Determined  Growth  Condition  (SDGC).  But  the  SDGC  may  not  hold  for  some 
infinite-dimensional  LTI  systems;  indeed  the  evolution  can  grow  exponentially  even  though  the  infinitesimal 
generator  (the  A-operator)  has  spectrum  inside  the  left  half  of  the  complex  plane  [6]-[8].  In  the  present  work  we 
use  perturbation  analysis  to  find  simple  conditions  under  which  the  spatially  periodic  system  satisfies  the  SDGC 
and  is  exponentially  stable. 

Our  presentation  is  organized  as  follows.  Section  II  outlines  the  main  results  of  the  paper.  Section  III  briefly 
reviews  the  frequency  representation  of  spatially  periodic  operators.  Section  IV  introduces  the  problem  setup. 
Section  V  discusses  the  analytic  perturbation  of  the  Tf^-norm  and  Section  VI  provides  related  illustrative  examples. 
Section  VII  studies  conditions  under  which  a  spatially  periodic  system  is  exponentially  stable.  Many  proofs  and 
technical  details  have  been  relegated  to  the  Appendix  to  improve  readability. 

Notation 

We  use  fc  G  M  to  characterize  the  spatial-frequency  variable,  also  known  as  the  wave-number.  S(r)  is  the 
spectrum  of  the  operator  T,  and  Sp(T)  its  point  spectrum,  p{T)  its  resolvent  set,  and  R{(,T)  its  resolvent 
operator  —  T)~^.  C“  denotes  all  complex  numbers  with  real  part  less  than  zero,  and  j  :=  \/—l.  denotes 
the  complex-conjugate  transpose,  and  also  the  adjoint  of  a  linear  operator.  6  is  the  closure  of  the  set  6  C  C.  We 
will  use  the  same  notation  A  for  an  operator  and  it’s  kernel  representation.  The  spatio-temporal  function  u{t,  x) 
(operator  A)  is  denoted  by  u{t,  k)  (respectively  A)  after  the  application  of  a  Fourier  transform  on  the  spatial 
variable  x.  Where  there  is  no  chance  of  confusion,  we  use  the  same  notation  for  a  spatially  invariant  operator 
and  its  Fourier  symbol. 

Terminology 

Throughout  the  paper,  we  use  the  terms  spatial  ‘operators’  and  spatial  ‘systems’.  By  the  former,  we  mean 
a  purely  spatial  system  with  no  temporal  dynamics  (i.e.  a  memoryless  operator  that  acts  on  a  spatial  function 
and  yields  a  spatial  function),  whereas  the  latter  refers  to  a  spatio-temporal  system  (a  system  with  an  internal 
state  evolves  on  some  spatial  domain,  i.e.,  for  every  time  t  the  state  is  a  function  on  a  spatial  domain).  When 
spatially  periodic  feedback  operators  are  small  in  norm,  we  will  use  the  phrases  periodic  feedback  and  periodic 
perturbation  interchangeably.  We  use  the  term  pure  point  spectrum  (or  equivalently  discrete  spectrum)  to  mean 
that  the  spectrum  of  an  operator  consists  entirely  of  isolated  eigenvalues  [9].  By  a  scalar  system,  we  mean  that 
the  Euclidean  dimension  of  the  state  is  equal  to  one. 

II.  Main  Results 

We  consider  systems  described  by  linear,  time-invariant,  integro  partial  differential  equations  defined  on  an 
unbounded  one  dimensional  domain.  We  use  a  standard  state-space  representation  of  the  form 

[dtf]{t,x)  =  [Af]{t,x)  -f  [Bu]{t,x), 
y{t,x)  =  [C'tp]{t,x), 


(1) 
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where  t  €  [0,  oo)  and  x  G  M.,  ip,  u,  y  are  spatio-temporal  functions,  and  A,  B,  C  are  spatial  integro-differential 
operators  with  coefficients  that  are  periodic  functions  with  a  common  period  X.  We  refer  to  such  systems  as 
spatially  periodic. 

It  is  often  physically  meaningful  to  regard  the  spatially  periodic  operators  as  additive  or  multiplicative  per¬ 
turbations  of  spatially  invariant  ones  [and  by  spatially  invariant  we  mean  integro-differential  operators  with 
constant  coefficients].  For  example,  the  generator  in  (1)  can  often  be  decomposed  as  A  =  A°  +  eE,  where 
A°  is  a  spatially  invariant  operator  and  E  is  an  operator  that  includes  multiplication  by  periodic  functions.  In 
some  control  applications,  the  operator  E  is  something  to  be  “designed”.  Therefore  it  is  desirable  to  have  easily 
verifiable  conditions  for  stability  and  norms  of  such  systems.  This  would  then  allow  for  the  selection  of  the  spatial 
period  and  amplitude  of  E  to  achieve  the  desired  behavior.  The  perturbation  analysis  we  present,  though  limited 
to  small  values  of  e,  provides  useful  results  for  selecting  candidate  “periods”  for  E. 

Our  analysis  and  results  are  derived  using  a  special  frequency  representation.  We  show  that  the  spatial  periodicity 
of  the  operators  A,  B  and  C  implies  that  (1)  can  be  rewritten  as 

[dtipe](t)  =  [A0ipe]{t)  +  [Beue]{t), 

ye{t)  =  [Ceip9]{t),  (2) 

where  9  G  [0,27r/2f);  for  every  value  of  9,  ipg,  ug,  yg  are  bi-infintie  vectors,  and  Ag,  Bg,  Cg  are  bi-infinite 
matrices.  The  systems  (2)  and  (1)  are  related  through  a  unitary  transformation,  and  in  particular  quadratic  forms 
and  norms  are  preserved  by  this  transformation.  Consequently,  stability  and  quadratic  norm  properties  of  (2)  and 
(1)  are  equivalent.  With  this  transformation,  the  analysis  of  the  original  system  (1)  is  reduced  to  that  of  the  family 
of  systems  (2)  that  are  decoupled  in  the  parameter  9.  In  particular,  perturbation  analysis  for  (2)  is  easier  and  less 
technical  than  that  for  the  original  system  (1). 

To  make  for  easier  reading  we  first  present  the  results  on  perturbation  analysis  of  the  Tf^-norm,  and  then  deal 
with  the  issue  of  stability. 

For  a  large  class  of  infinite-dimensional  systems,  computing  the  Tf^-norm  involves  solving  an  operator  algebraic 
Lyapunov  equation 

AP  +  PA*  =  -BB*. 

In  general  this  is  a  difficult  task  that  must  be  done  using  appropriate  discretization  techniques.  However,  in  the 
case  when  A  and  B  are  spatially  periodic  operators,  then  so  is  the  solution  P.  Thus  the  frequency  representation 
implies  that  this  operator  Lypunov  equation  is  equivalent  to  the  decoupled  family  of  matrix  Lyapunov  equations 

AgVg  +  V0A*e  =  -BgB*e,  (3) 

where  Ag,Bg  and  Vg  are  the  bi-infinite  matrix  representations  of  A,  B  and  P.  Once  Vg  is  found,  the  Tf^-norm 
of  the  system  can  be  computed  from  [5] 

iv&ce[CeVg C*g]  d9,  O  =  27r/X. 

Solving  (3)  is  still  a  difficult  problem  in  general  since  it  involves  bi-infinite  matrices.  We  use  perturbation 
analysis  as  follows:  the  generator  is  expressed  as  Ag  =  Ag  +  eSg  where  Ag  and  Eg  correspond  to  the  spatially 
invariant  and  spatially  periodic  components  respectively.  It  follows  that  the  solution  Vg  is  analytic  in  e  and  the 
terms  of  its  power  series  expansion  (denoted  by  Vg''^)  satisfy  a  sequence  of  forward  coupled  Lyapunov  equations. 
Furthermore,  the  terms  v'g'^  are  banded  matrices,  with  the  number  of  bands  increasing  with  the  index  i.  These 
Lyapunov  equations  can  then  be  solved  recursively  for  i  =  0,1,2,....  Formulae  for  these  representations  and 
the  corresponding  sequence  of  Lyapunov  equations  are  given  in  Section  V.  In  some  examples  that  we  present 
in  Section  VI,  these  formulae  lead  to  simple  “resonance”  conditions  for  stabilization  or  destabilization  of  PDFs 
using  spatially  periodic  feedback. 

The  second  set  of  results  concern  the  problem  of  stability.  As  mentioned  in  the  Introduction,  when  A  is 
an  infinite-dimensional  operator  it  is  possible  that  its  spectrum  S(A)  lies  inside  C“,  and  yet  ||e^‘||  grows 
exponentially  [6]-[8].  In  such  cases  it  is  said  that  the  spectrum-determined  growth  condition  is  not  satisfied  [8]. 
Yet  there  exists  a  wide  range  of  infinite-dimensional  systems  for  which  the  spectrum-determined  growth  condition 
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Fig.  1.  Left:  Pictorial  representation  of  A{k,  k)  ■*/>(«;)  dK.  Right:  The  frequency  kernel  A  of  a  spatially  periodic  operator  A. 


is  satisfied.  These  include  (but  are  not  limited  to)  systems  for  which  the  A-operator  is  sectorial  (also  known  as 
an  operator  which  generates  a  holomorphic  or  analytic  semigroup)  [9]-[ll]  or  is  a  Reisz-spectral  operator  [2].  In 
this  paper  we  focus  on  sectorial  operators. 

Thus  to  establish  exponential  stability  of  a  system,  one  possibility  would  be  to  show  simultaneously  that 

(i)  A  is  sectorial, 

(ii)  I](A)  lies  in  C“. 

But  this  still  does  not  make  the  problem  trivial.  In  fact  proving  that  an  infinite-dimensional  operator  is  sectorial, 
and  then  finding  its  spectrum,  can  in  general  be  extremely  difficult. 

Once  again  we  use  perturbation  methods  to  show  (i)  and  (ii).  We  consider  A  to  have  the  form  A  =  A°  +  cE 
where  A°  is  a  spatially  invariant  operator,  E  is  a  spatially  periodic  operator,  and  e  is  a  small  complex  scalar. 
Using  the  bi-infinite  matrix  representation,  we  first  find  conditions  on  the  spatially  invariant  operator  A°  such 
that  (i)  and  (ii)  are  satisfied.  We  then  show  that  (i)  and  (ii)  will  remain  satisfied  if  e  is  small  enough  and  if  the 
spatially  periodic  operator  E  is  ‘weaker’  than  A°  (in  the  sense  that  E  is  relatively  bounded  with  respect  to  A°). 
The  utility  of  this  approach  is  that  (i)  and  (ii)  are  much  easier  to  check  for  a  spatially  invariant  operator  than  they 
are  for  a  spatially  periodic  one. 

III.  Frequency  Representation  oe  Periodic  Operators 
In  this  section  we  briefly  discuss  the  frequency  domain  representation  of  spatially  periodic  operators.  We  then 
show  how  this  representation  can  be  used  to  convert  a  spatially  periodic  system  into  a  family  of  matrix-valued 
LTI  systems.  For  a  detailed  account  the  reader  is  referred  to  [5]  and  [12]. 

Let  ip{k)  and  <j>{k)  denote  the  Fourier  transforms  of  two  spatial  functions  'ijj{x)  and  <j>{x)  respectively.  If  and 
<j)  are  related  by  a  linear  operator,  (j)  =  Aip,  then  so  are  ip  and  p  and  we  have 

p{x)  =  f  A{x,x)'P{x)  dx  ^  (p{k)  =  f  A{k,K)ip{n)  dn  (4) 

where  A  and  A  are  kernel  functions  in  the  spatial  and  Fourier  domains,  respectively;  see  Figure  1  (left).  It  is 
shown  in  [5]  [12]  that  the  most  general  spatially  periodic  operator  A  with  spatial  period  X  =  27r/U  can  be 
represented  in  the  Fourier  domain  as  an  operator  with  a  kernel  function  of  the  form 

A{k,K)  =  Ai{k)  S{k  —  K  —  fll).  (5) 

/ez 

Thus  the  kernel  function  corresponding  to  A  is  composed  of  parallel  and  equally-spaced  ‘impulse  sheets’  which 
can  be  visualized  in  Figure  1  (right).  [5]  further  describes  how  the  particular  structure  (5)  of  A  allows  (4)  to  be 
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given  a  matrix  representation 


Ao(e-r!)  A_i(e-ct)  A_2(e-n) 

Aim  Aom  A_i(e) 

A2(e+n)  Ai(»  +  ri)  A(](e  +  S2) 

for  which  we  adopt  the  notation 

(jig  =  Ag  'ipg. 

In  other  words,  a  general  spatially  periodic  operator  A  can  be  described  by  a  family  of  (bi-infinite)  matrices  Ag 
parameterized  by  a  variable  6. 

Remark  1:  We  emphasize  that  the  Ai{-)  in  (5)  and  (6)  can  be  matrices.  Thus,  in  general,  Ag  has  a  “block” 
structure.  But  throughout  this  paper  and  for  the  sake  of  simplicity,  we  choose  not  to  explicitly  refer  to  this  block 
structure,  even  though  all  our  results  hold  for  matrix-valued  • ).  In  the  same  light,  we  do  not  refer  directly  to 
the  Euclidean  dimension  of  the  vectors  ■ )  and  </>(•).  ■ 

Spatially  invariant  [13]  and  spatially  periodic  pure  multiplication  operators  constitute  special  subclasses  of 
spatially  periodic  operators.  In  the  framework  established  above,  Ag  is  diagonal  for  spatially  invariant  operators, 
and  Toeplitz  for  periodic  pure  multiplication  operators. 

Example  1:  A  =  dx  and  F{x)  =  cos{Elx)  have  the  following  representations 


Ag  =  jd+jQ.n  ,  T 


for  every  9  G  [0,17),  respectively.  Notice  that  since  A  is  spatially  invariant,  it  is  fully  described  by  its  Fourier 
symbol  A[){k)  =  jk,  fc  G  M.  And  it  is  the  samples  of  Ao(  • )  that  make  up  the  diagonal  of  Ag  for  every  9.  We 
have  dropped  the  9  subscript  in  T,  as  it  is  independent  of  this  variable.  ■ 

Remark  2:  It  is  possible  to  define  a  unitary  operator  ^g  [5]  such  that  ipg  =  ^g  ij},  (pg  =  ^g  cp,  and  thus 
Ag  =  ^g  A^^ .  is  equivalent  to  the  frequency  domain  lifting  operation  of  [14]  and  [15]  (see  also  [16]).  By 
the  unitary  property  of  the  lifting  operator  it  follows  that 

p  p^ 

/  ivdJce[Ai{k)  Alfk)]  dk  =  /  tr ace[ Ag  Ag]  d9  =  /  |j^e|lHs 

/gZ  JR  Jo  Jo 

with  |jT||^g  :=  trace[TT*]  being  the  square  of  the  Hilbert-Schmidt  norm^  of  T.  ■ 

Finally,  given  a  spatially  periodic  system  in  state-space  form  (1)  with  spatially  periodic  A,  B,  and  C,  one  can 
replace  each  of  these  operators  with  its  bi-infinite  matrix  representation  to  obtain  the  family  LTI  systems  (2). 

Example  2:  Consider  the  spatially  periodic  heat  equation  on  x  G  M 

dtf{t,x)  =  {dl  —  a  cos(nx))  ip{t,  x)  +  u{t,  x), 

y{t,x)  =  'ip{t,x),  (8) 

with  real  a  ^  0  and  12  >  0.^  Clearly  A  =  +  a  cos{idx)  with  domain 

^  =  {(p  £  L^CR)  I  (p,  ^  absolutely  continuous,  G  L^(IR)}, 
dx  dx^ 

^The  Hilbert-Schmidt  norm  of  an  operator  is  a  generalization  of  the  Frobenius  norm  of  finite-dimensional  matrices  ||A||p  = 

E™.n  kmnp  =  trace[AA*]. 

^By  dtipit^x)  and  we  mean  the  spatio-temporal  functions  dt'ijj  and  d^ip  evaluated  at  the  point  {t,x). 


if(O-Q) 

^(0)  ,  OG[o,n),  (6) 

4>{o  +  n) 
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B  =  C  =  6{x)  are  the  identity  convolution  operator  on  L^(IR).  Rewriting  the  system  in  its  matrix  representation 


dttpeit)  =  Aeipeit)  +  B0ug{t), 
yg{t)  =  Cg-ipeit), 

with 


r  n 

r  -1 

Ag  = 

—  (0  -f  17n)^ 

—  a 

■■■  0  5 

1  0  ■■■ 

,  Bg  =Cg  = 

1 

GeH  :=  Cg  {jujJ- Ae)-^Bg 


-1 


CK 

2 

f  ju;  +  (0+nn)2  I 

a 

2 


(9) 


(10) 


(11) 


Notice  that  (9)-(10)  is  now  fully  decoupled  in  the  variable  6.  In  other  words,  (8)  is  equivalent  to  the  family  of 
state-space  representations  (9)-(10),  and  transfer  function  (11),  all  parameterized  by  the  variable  9  €  [0,17).  ■ 


IV.  Problem  Setup 

Let  us  now  consider  a  system  of  the  form 

dttp{t,x)  =  Ail;{t,x)  +  Bu{t,x) 

=  {A°  +  B°  eF  C°)ilj{t,x)  +  Bu{t,x),  (12) 

y{t,x)  =  C'tf;{t,x), 

where  t  G  [0,  oo)  and  a;  S  M  with  the  following  assumptions.  The  (possibly  unbounded)  operators  A°,  B°,  C°  are 
spatially  invariant,  and  the  bounded  operators  B,  C  are  spatially  periodic.  F{x)  =  2Lcos{^lx)  with  L  a  constant 
matrix,  and  e  is  a  complex  scalar.  A°,  B°,  C°  and  E  :=  B°  F  C°  are  all  defined  on  a  dense  domain  Si  C 
and  are  closed,  u,  y,  and  i/;  are  the  spatio-temporal  input,  output,  and  state  of  the  system,  respectively.  We  will 
refer  to  A  as  the  infinitesimal  generator  of  the  system. 

Comment  on  Notation:  To  avoid  clutter,  we  henceforth  drop  the  “ '  ”  on  the  Fourier  symbol  of  operators  and 
frequency  domain  functions.  For  example,  we  use  Ao(  • )  [instead  of  Ao(  • )]  to  represent  the  Fourier  symbol  of 
the  spatially  invariant  operator  A°.  m 

As  shown  in  detail  in  [5]  and  also  briefly  in  the  previous  section,  the  system  can  be  represented  in  the  (spatial) 
Fourier  domain  by  the  family  of  systems 

dtlpeit)  =  Agfig{t)  +  BgUg{t) 

=  {A°g  +  eB°gFC°g)  Mt)  +  Bg  Ug{t)  (13) 

yg{t)  =  Cglpgit), 

parameterized  by  0  S  [0,17).  Here  Bg  and  Cg  have  the  general  form  of  the  operator  in  (6)  [i.e.  can  possess  any 
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Fig.  2.  has  a  spatially  invariant  infinitesimal  generator  A°.  The  LFT  of  and  the  spatially  periodic  multiplication  operator  eF  yields 
a  system  which  has  a  spatially  periodic  infinitesimal  generator  A. 


number  of  nonzero  sbudiagonals],  T  has  the  form  given  in  Example  1  with  \  replaced  by  L,  and 


^o(^n) 

,  B°o  = 

BAG  A 

po  — 

5  '^9  — 

CAOA 

^2  = 


with  :=  0  +  rin,  n  G  Z,  and 

Ai(-)  :=  B°{-)LC°{  -  -n), 


So  ■.=  BlTa  = 


0  A_i(0„) 
Ai(e„)  0 


(14) 


A_i(-)  :=  B°{-)LC°{  -  +n). 


(15) 


We  emphasize  that  the  convention  used  in  the  representation  of  So  in  (14)  is  the  same  as  that  used  in  (6);  for 
example  the  n*  row  of  is  {•  •  •  ,  0,  Ai(0„),  0,  A_i(0„),  0,  •  •  •  }. 

Remark  3:  We  note  that  taking  F{x)  to  be  a  pure  cosine  is  not  restrictive.  The  results  obtained  here  can 
be  easily  extended  to  problems  where  F{x)  is  any  periodic  function  with  absolutely  convergent  Fourier  series 
coefficients.  The  inclusion  of  higher  harmonics  of  St  in  F{x),  namely  functions  of  frequency  2St,  3St,  etc,  would 
not  reveal  new  interesting  phenomena  and  would  only  complicate  the  algebra.  ■ 

Remark  4:  The  system  (12)  can  be  considered  as  the  LFT  (linear  fractional  transformation  [17])  of  a  spatially 
periodic  system  G°  with  spatially  invariant  infinitesimal  generator  A°, 


B 

B° 

G°  = 

G 

0 

0 

.  c° 

0 

0 

and  the  (memoryless  and  bounded)  spatially  periodic  pure  multiplication  operator  eF{x)  =  e2Lcos{Stx),  see 
Figure  2.  ■ 


Stability  Analysis  and  Sectorial  Operators 

It  is  shown  in  [5]  that  for  a  general  spatially  periodic  operator  A  we  have 

S(^)  =  U  Z{Ao).  (16) 

ee[o.a) 

In  the  case  where  A  is  spatially  invariant  [and  thus  Ao  =  diag  {•  •  •  ,  Ao(0„),  •••}],  (16)  further  simplifies  to 

S(^)  =  USp(Ao(fc)).  (17) 

keE 

Example  3:  Let  A  =  —{d^  +  Then  Ao{k)  =  —{k^  —  see  Figure  3  (left).  On  the  other  hand,  since 
Aq{  ■ )  is  scalar,  Sp(Ao(fc))  =  AQ{k)  for  every  k.  It  is  easy  to  see  that  ylo(  • )  takes  every  real  negative  value  and 
thus  from  (17)  A  has  the  continuous  spectrum  S(A)  =  (— oo,0],  see  Figure  3  (center).  ■ 

Remark  5:  When  A  is  spatially  invariant,  a  helpful  way  to  think  about  S(^)  in  terms  of  its  symbol  Aq  is 
suggested  by  the  previous  example.  First  plot  Sp(ylo(  • ))  in  the  ‘complex-plane  x  spatial-frequency’  space,  as  in 
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Fig.  3.  Left:  Representation  of  the  symbol  Ao(  • )  of  Example  3  in  ‘complex-plane  x  spatial-frequency’  space.  Center:  S(A)  in  the  complex 
plane.  Right:  For  spatially  invariant  A,  the  (diagonal)  elements  of  J^q  are  samples  of  the  Fourier  symbol  Ao(  • ). 


Figure  3  (left)  of  Example  3.  Then  the  orthogonal  projection  onto  the  complex  plane  of  this  plot  would  give  S(A), 
as  in  Figure  3  (center).  This  can  be  considered  as  a  geometric  interpretation  of  (17).  In  Example  3,  since  ^o(  • ) 
is  real  scalar  and  takes  only  negative  values,  this  projection  yields  only  the  negative  real  axis.  But  in  general 
if  ^o(')  G  this  projection  would  lead  to  q  curves  in  the  complex  plane.  Notice  also  that  in  this  setting, 

Ti{Ag)  is  the  projection  onto  the  complex  plane  of  samples  of  Sp(Ao(  • ))  taken  at  fc  =  =  0  +  fin,  n  G  Z, 

in  the  ‘complex-plane  x  spatial-frequency’  space.  This  can  be  considered  as  a  geometric  interpretation  of  (16). 
Figure  3  (right)  shows  these  samples  in  the  ‘complex-plane  x  spatial-frequency’  space  for  a  scalar  A.  ■ 

We  next  introduce  a  special  subclass  of  holomorphic  (or  analytic)  semigroups.  The  reader  is  referred  to  [9]- 
[11]  for  a  detailed  discussion.  Suppose  A  is  densely  defined,  p{A)  contains  a  sector  of  the  complex  plane 
I  arg(z  —  a)  I  <  ^  -f  7,  7  >  0,  a  G  M,  and  there  exits  some  M  >  0  such  that 

\\{zI-A)-^\\  <  for  |arg(z-a)|  <  ^  +  7-  (18) 

\z  —  a\  I 

Then  A  generates  a  holomorphic  semigroup  and  we  write  A  G  J^(7,a,  M)  [11]  [9].  We  say  that  A  is  sectorial 
if  A  belongs  to  some  a,  M). 

Finally,  a  semigroup  is  called  exponentially  stable  if  there  exist  positive  constants  M  and  g  such  that  [2] 

\\e^^\\  <  for  t  >  0. 

Theorem  1:  Fet  A  be  sectorial.  Then  if  S(Gl)  C  C“,  A  generates  an  exponentially  stable  semigroup. 

Proof:  If  A  is  sectorial  it  defines  a  holomorphic  semigroup  and  thus  is  differentiable  for  t  >  0  [10].  Then 
[8]  shows  that  this  is  sufficient  for  the  spectrum-determined  growth  condition  to  hold.  In  particular,  if  S(g1)  C  C“, 
A  generates  an  exponentially  stable  semigroup.  ■ 


Tif-Norm  of  Spatially  Periodic  Systems 

We  define  the  Tf^-norm  of  an  exponentially  stable  spatially  periodic  system  G  as 


=  I 


tvdiCe[Qe{t)  Gg  (f)]  dt  dO, 


0  ^0 


where  Gg{t)  =  Cge^^^Bg  is  the  impulse  response  of  the  system  (13).  The  intuition  for  this  definition  and  the 
proof  of  the  following  theorem  can  be  found  in  [5]. 

Theorem  2:  Consider  the  exponentially  stable  spatially  periodic  FTI  system  G,  with  spatial  period  X  =  2tt/GI 
and  state-space  realization  (12)-(13).  We  have 

\\G\\n-  =  ^  !  t^txce[CgVgCl]de  =  ^  [  tra.ce[B;  QgBg]de, 

27r  Jo  Jo 

where  Vg  and  Qg  are  the  solutions  of  the  (0-parameterized)  algebraic  Fyapunov  equations 


Ag  Qg  +  Qg  Ag  —  —Cg  Cg . 


AgVg  +  VgA*g  =  -BgB*g, 
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V.  Perturbation  Analysis  of  the  TY^-Norm 

The  difficulty  in  calculating  the  Tf^-norm  using  Theorem  2,  is  that  unless  Ag,  Bg  and  Cg  are  diagonal  (i.e.,  G 
is  a  spatially  invariant  system),  the  operators  Vg  and  Qg  are  “full”,  meaning  that  they  possess  all  of  their  (infinite 
number  of)  subdiagonals.  This  makes  the  computation  of  the  Tf^-norm  numerically  intensive.  Namely,  one  has 
to  solve  an  infinite-dimensional  algebraic  Lyapunov  equation  to  find  the  operator  Vg  (or  Qg)  for  every  value  of 
9  €  [0,  n).  In  this  section  we  will  see  how  one  can  use  analytic  perturbation  techniques  to  compute  the  Tf^-norm 
in  a  simple  and  numerically  efficient  way,  and  without  having  to  explicitly  find  the  full  Vg  and  Qg  operators. 
Such  a  perturbation  analysis  is  very  useful  in  predicting  general  trends  and  extracting  valuable  information  about 
the  Tf^-norm. 

Let  us  now  consider  the  general  setup  of  (12),  where  we  take  e  is  a  small  real  scalar.  We  also  assume  that  both 
A°  and  A  =  A°  +  eE  define  exponentially  stable  strongly  continuous  semigroups  (also  known  as  Co-semigroups) 
on  L^(IR)  [2],  and  that  B  and  C  are  spatially  invariant  operators.  We  are  interested  in  the  changes  in  the  Tf^-norm 
of  this  system  for  small  magnitudes  of  e  and  different  values  of  the  frequency  fl. 

Let  us  define 

Vg{e)  :=  iP^°)  +  eipW+e2iPf)+...  , 

with  Vg{e)  =  Vg{e).  Notice  that  this  implies  V^'’*  =  V^^  for  all  m  =  0, 1, 2,  •  •  • .  Our  aim  is  to  find  7^^™^  by 
solving  the  Lyapunov  equation 

Ag{e)  Vg{e)  +  Vg{e)A*g{e)  =  -Bg  (19) 

{A°g  +  eSg)  {Vf  ^  +  eV^P  +  vP  +  •••)  + 

{vP  +  tVP  +  e"  vP  +  •  •  • )  (-4^  +  =  -Bg  B*g, 

and  compute  the  Tf^-norm  of  the  system  using  Theorem  2 

1 

\\G\\u^  =  TT  /  A^ce[CgVg{e)C*g\de. 

It  is  easy  to  see  from  (20)  that 

A°gVP  +  vP  AT  =  -BgBT 
A°gVP  +  vP  AT  =  -pgVP  +  vPs*g), 

A^gvP  +  vP  AT  =  -pevP  +  vPs;), 

Now  since  A°g  and  Bg  Bg  are  diagonal  in  (21),  so  is  Vp\  In  (22),  the  right  hand  side  operator  has  the  structure  of 
being  nonzero  only  on  the  first  upper  and  lower  subdiagonals,  and  hence  vP  inherits  the  same  structure  (since 

(n\ 

A°g  is  diagonal).  In  the  same  manner,  one  can  show  that  Vg  '  is  only  nonzero  on  the  main  diagonal  and  the  second 
upper  and  lower  subdiagonals,  and  so  on  for  other  Vp\  We  have 


vP  =  Po{eA  ,  vP  =  '■  °  ,  vP 

PiidA  0 
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where  the  n*  row  of  is  equal  to  {•••  ,  0,  Pi(0„),  0,  Pf  (0„),  0,  •  •  •  },  and  the  n*  row  of  is  equal  to 

{•  •  • ,  0,  P2(0„),  0,  go(0„),  0,  P2*(0„),  0,  •  •  • }. 

Remark  6:  It  is  important  to  realize  that,  not  only  is  not  a  “full”  operator,  it  has  at  most  m  nonzero 

upper  and  lower  subdiagonals.  Also,  all  for  odd  m  have  zero  diagonal  and  are  thus  trace-free  operators.  ■ 
Remark  7:  Note  that  even  though  Ag  =  Ag  +  eSg  has  only  one  nonzero  subdiagonal  [see  (14)],  Pe(e)  = 
+  eVg^'^  +  ■  ■  ■  possesses  all  of  its  subdiagonals.  This  is  precisely  the  reason  why  direct  calculation  of  Vg  in 
Theorem  2  is  computationally  difficult.  ■ 

Now  returning  to  (21)-(23),  Vg^\  Vg^'^  and  V^g^  are  found  by  equating,  element  by  element,  the  bi-infinite 
matrices  on  both  sides  of  these  equations.  For  example,  (21)  leads  to 

Aq{6  +  fin)  Po{0  +  fin)  +  Po{0  +  fin)  Aq{0  +  fin)  =  —  B (9  +  fin)  B*  (9  +  fin) 

for  every  n  G  Z,  and  9  G  [0,12).  But  notice  that  as  n  assumes  values  over  all  integers  and  9  changes  in  [0,12), 
k  =  9  +  fin  takes  all  real  values,  and  one  can  rewrite  the  above  equation  as 

Ao{k)Po{k)  +  Po{k)A*{k)  =  -B{k)B*{k), 

for  all  fc  G  M.  Applying  the  same  procedure  to  (22)-(23),  one  arrives  at 


Ao{k)  Po{k)  +  Po{k)  A*{k)  =  -B{k)  B*{k),  (24) 

Ao{k)  Pi{k)  +  Pi{k)  A;{k  -  fl)  =  -(^Ai{k)Po{k-fl)  +  Po{k)A*_^{k-fl)y  (25) 

Ao{k)Qo{k)  +  Qo{k)A*{k)  =  -(^A,{k)P*{k)  +  Pi(fc)A*(fc)  +  A_,{k)  P,{k  +  fl) 

+  P*{k  +  fl)A*_,{k)y  (26) 

and  so  on  for  all  nonzero  diagonals  of  Vg  ,  m  =  3, 4,  •  •  • . 

Remark  8:  Notice  that  from  the  above  equations,  one  first  finds  Po(  • )  from  (24),  then  Pi(  • )  from  (25),  and 
so  on.  In  other  words,  computing  the  subdiagonals  of  Vg  becomes  “decoupled”  in  one  direction.  This  decoupling 
would  not  have  been  possible  had  we  not  employed  a  perturbation  approach  and  had  attempted  to  solve  (19) 
directly.  ■ 

Returning  to  the  calculation  of  the  Tf^-norm,  let  us  first  separate  the  diagonal  part  of  by  rewriting  it  as 


P^2) 


where 


vf 


Qo{dn 


and  Pp^  contains  the  rest  of  V^\  Clearly  trace [Cg  Pp^  =  0.  Also,  recall  that 

tracepe  pp^+i)  =  0,  m  =  0, 1, 2,  •  •  •  . 
Now  one  can  write  the  following 


(27) 


IIGII?,.  = 


2it 


tracep^  V0{e)  C^]  dO 


1 

—  /  trace[Ce  Pp^  Cg  vf^  C*g]  d9  +  0(ep 

Jo 

[  tra.ce[CgVj°^C*g  +  t^CgVfc*g]d9  +  0(ep, 

Jo 


where  the  absence  of  odd  powers  of  e  is  due  to  (27),  and  the  last  equation  follows  from  the  fact  that  trace  [Cg  Vg  Cg] 
0.  Next,  using  the  unitaryness  of  the  lifting  transform  we  have 

pVL  poo 

/  trace[CePp^Q]d6»  =  /  trace [C'(fc)  Po(fc)  C'*(fc)]  dfc, 

J  0  _  J  —  oo 

poo 

/  iv&ce[CgVg  C*g]d9  =  /  ivace[C{k)Qo{k)C*{k)]dk, 

Jo  J  —  oo 
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and  we  arrive  at  the  final  result 
1  1'°° 

\\G\\^2  =  —  tiace[C{k)Po{k)C*{k)  +  e'^C{k)Qo{k)C*{k)]dk  +  0{e^).  (28) 

27r 

Remark  9:  Let  us  stress  again  the  important  advantage  of  the  above  perturbation  analysis  in  comparison  to 
the  method  of  Theorem  2.  In  the  direct  method  of  computing  the  Tf^-norm  proposed  in  Theorem  2,  one  has  to 
solve  a  family  of  infinite-dimensional  (operator)  algebraic  Lyapunov  equations,  whereas  the  perturbation  method 
reduces  the  Tf^-norm  computation  to  that  of  solving  a  family  of  finite-dimensional  (matrix)  Lyapunov/Sylvester 
equations  (24)-(26).  ■ 

In  summary,  we  have  the  following. 

Main  Result  of  Perturbation  Analysis  of  Ti.^ -Norm:  Consider  the  exponentially  stable  spatially  periodic  LTI 
system  G,  with  spatial  period  X  =  27r/n  and  state-space  realization  (12).  Then  for  small  values  of  |e|  the  Ti.^- 
norm  of  the  system  (12)  can  be  computed  from  (28),  where  Po{-)  and  Qo(  ■)  are  solutions  of  the  family  of 
finite-dimensional  Lyapunov/Sylvester  equations  described  by  (24)-(26).  ■ 


VI.  Examples 

As  an  application  of  the  perturbation  results  of  the  previous  section,  we  first  investigate  the  occurrence  of 
‘parametric  resonance’  for  a  class  of  spatially  periodic  systems.  Parametric  resonance  occurs  when  a  specific 
frequency  12^8  of  the  periodic  perturbation  resonates  with  some  ‘natural  frequency’  of  the  unperturbed  system, 
leading  to  a  local  (in  fl)  change  in  system  behavior  [3].  In  the  systems  we  consider  in  this  section,  this  change 
in  behavior  is  captured  by  the  value  of  the  Tf^-norm. 

Example  4:  Let  us  consider  the  periodic  PDE 

dtf  =  —  c  fit  -\-  e  cos{flx)  4-  u 

y  =  f,  (29) 


with  0  K  G'R  and  c  >  0.  Comparing  (29)  and  (12)  we  have 

AQ{k)  = -{k"^  -  xY  -  c,  B°{k)  =  l,  G°{k)  =  l,  B{k)  =  l,  C{k)  =  l, 
Eor  this  system,  the  functions  Po{k)  and  Qo{k)  of  the  previous  section  simplify  to^ 


Po{k) 

Qo{k) 


-1 

2Ao(fc) 

1  /  -1  -1  \ 

(Ao(fc))^  '2Ao(fc  —  n)  2Ao(fc-l-r2)/ 

4{Po{k))\Po{k  -n)  +  Po{k  +  n)), 


and  it  is  our  aim  to  find  the  Tf^-norm 

1  1"°° 

\\Gfn-  =  —J  (Poik)  +  e^Qo{k))dk  +  0{e^), 


(30) 

(31) 


(32) 


for  different  values  of  the  parameter  12  >  0.  More  specifically,  we  are  interested  in  the  values  of  12  for  which  the 
Tf^-norm  is  maximized. 

Erom  (30),  Po{k)  =  |  (fc2_^2)2_|_e-  The  first  plot  of  Eigure  4  shows  Po{  ■ ),  while  the  second  shows  Po(  •  ~ 
and  Po{-  4-12)  (dashed),  for  a  given  value  of  12  0.  As  12  is  increased,  Po(  •  —  12)  slides  to  the  right  and 

Pq{  •  4- 12)  to  the  left.  Erom  (31)  it  is  clear  that  to  find  Qo{  ■ )  for  any  12,  one  would  sum  the  two  functions  in 
the  second  plot  and  multiply  the  result  by  the  square  of  the  first  plot.  The  interesting  question  now  is,  for  what 
value(s)  of  12  S  (0,  oo)  would  the  ^-norm  in  (32)  be  maximized. 


^To  find  Qo(k)  one  needs  to  first  find  Pi(fc),  but  we  have  omitted  the  details  for  brevity. 
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Fig.  4.  Above:  Plot  of  Po(  ■  )•  Below:  Plot  of  Po(  ■  —  and  Po{-  +  (dashed). 


6 


Fig.  5.  Left:  Plots  of  Example  4  for  =  1  and  c  =  0.1.  Notice  that  the  first  graph  is  plotted  against  k  and  the  second  against  fl.  Right: 
The  plot  of  H^-norm  of  the  same  example,  but  calculated  by  taking  large  truncations  of  the  Ae,  ^9,  Cq  matrices  and  using  Theorem  2. 


One  can  easily  see  that  as  O  ^  0,  the  peaks  of  Po(  •  ~  Pq{  ■  +  O)  merge  toward  those  of  (Po(  •  ))^. 

thus  Qo{k)  dk  grows  and  hence  ||G||^2  grows.'*  This  is  not  surprizing;  as  O  ^  0,  the  perturbation  is  tending 
toward  a  constant  function  F{x)  =  cos(Ox)  ^  1.  This  results  in  shifting  the  whole  spectrum  of  A°  toward  the 
right-half  of  the  complex  plane,  thus  increasing  the  Ti^-norm  of  the  perturbed  system. 

But  we  are  more  interested  in  frequencies  0^0  that  lead  to  a  local  (in  O)  increase  in  the  Ti^-norm.  Now 
notice  that  a  different  alignment  of  the  peaks  can  also  occur,  which  leads  to  another  local  maximum  of  the  T-p- 
norm  as  a  function  of  fl.  This  happens  when  the  peak  of  Po(  •  —  fl)  at  fc  =  —k  +  becomes  aligned  with  the 
peak  of  (Po(  • ))  at  A:  =  K,  and,  simultaneously,  the  peak  of  Po(  •  -1-17)  at  k  =  x  —Q.  becomes  aligned  with  the 
peak  of  (Po(  •  ))^  at  k  =  —K.  Clearly  this  occurs  when 

— t”  Hjgg  —  —  2  . 

This  result  agrees  exactly  with  that  obtained  in  [12],  where  in  the  analysis  of  the  same  problem  it  is  shown  that  the 
part  of  the  spectrum  of  A  closest  to  the  imaginary  axis  ‘resonates’  with  perturbations  whose  frequency  satisfies 
the  relation  17  =  2  x:. 

Consider  (29)  with  x  =  1  and  c  =  0.1.  For  this  system  Po{k)  dk  «  4.74.  Figure  5  (left)  shows  the  graphs 
of  Po(  • )  (plotted  against  k)  and  Qo{k)  dk  (plotted  against  17).  The  peak  at  17  =  2  in  the  lower  plot  verifies 
the  result  17res  =  2  x  obtained  above. 

Figure  5  (right)  shows  the  Tf^-norm  of  this  system  computed  by  taking  large  enough  truncations  of  the  Ae, 
Bg,  Cg  matrices  and  then  applying  Theorem  2.  The  figure  shows  that  the  trends  were  indeed  correctly  predicted 
by  the  perturbation  analysis;  the  peaks  at  17  =  0,2  correspond  to  those  of  Qo(k)  dk. 

‘^Remember  that  Po{k)  is  independent  of  f2,  and  thus  Po{k)  dk  remains  constant  for  different  O. 
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Now  consider  (29)  with  x  =  \  and  c  =  0.1,  but  with  e  replaced  by  ej  (i.e.,  a  purely  imaginary  purturbation). 
Obviously  the  unperturbed  system  remains  the  same  as  before  and  hence  Po{k)  dk  «  4.74.  Figure  5  shows 
the  graph  of  Qo{k)  dk,  which  demonstrates  that  for  this  system  the  purely  imaginary  perturbation  reduces 
the  Tf^-norm  at  all  frequencies.  We  address  the  physical  interpretation  of  such  a  perturbation  in  the  Appendix.  ■ 

We  continue  with  more  examples  to  demonstrate  that  by  appropriately  choosing  the  frequency  of  the  perturba¬ 
tion,  one  can  decrease  or  (as  in  the  previous  example)  increase  the  Tf^-norm  of  the  unperturbed  system. 


Example  5:  Let  us  consider  the  periodic  PDF  [5] 

dtip  =  —{dl  +  x'^Y'ijj  —  c-ip  +  e  cos{flx)  dx  4’  +  u 

y  =  'ip. 


(33) 


Comparing  (33)  and  (12)  we  have 

Ao{k)  = -{k^  -  x^f  -  c,  B°{k)  =  l,  C°{k)=jk,  B{k)  =  l,  C{k)  =  l,  L  =  ^. 

The  difference  between  this  system  and  (29)  is  that  here  C°  is  a  spatial  derivative.  The  plot  of  Figure  7  demonstrates 
f^ao  Qo{k)  dk  for  =  1  and  c  =  0.1.  Notice  that  the  peak  at  =  2  remains  the  same  as  in  Figure  5,  but 
we  now  have  a  decrease  at  small  frequencies.  This  is  due  to  the  derivative  operator  C°  =  dx-  Finally,  notice  the 
agreement  of  the  perturbation  analysis  here  with  the  (non-perturbation)  calculations  for  the  same  system  in  the 
Third  Example  of  Section  VII  in  [5].  Our  perturbation  methods  correctly  predict  the  increase  at  O  =  2  and  the 
decrease  around  O  «  0.4  of  the  Tf^-norm.  ■ 

Example  6:  The  following  system  is  motivated  by  channel  flow  problems 


Ao{k)  = 


—  j^{k^  +  c)  0 

jk  -A{P  +  c) 


B°{k)  = 


1  0 
oil 


C°{k)  = 


1  0 
0  1  I 


B{k)  = 


1  0 
0  oj 


C{k)  = 


0  0] 
0  1  I 


We  perform  numerical  calculations  for  i?  =  6,  c  =  1.  We  have  ir&ce[C(k)PQ(k)C* (k,)]  dk  «  20.72.  Figure  8 
shows  that  the  Tf^-norm  can  be  decreased  by  the  the  application  of  periodic  perturbations  with  frequency  O  «  0.7. 
It  is  interseting  that  if  one  uses  the  locations  of  the  peaks  in  the  first  plot  of  Figure  8  to  find  x  =  0.75,  then 
from  the  peak  at  fires  =  1.6  in  the  second  plot  it  seems  that  the  relationship  fires  ~  2><'  =  1.5  still  holds  with  an 
acceptable  error  even  for  this  non-scalar  example.  ■ 
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Fig.  8.  Graphs  of  Example  6  for  i?  =  6  and  c  =  1. 


VII.  Stability  and  the  Spectrum-Determined  Growth  Condition 

In  the  literature  on  semigroups,  there  exist  examples  in  which  S(4l)  lies  entirely  inside  C“,  but  ||e"^*||  does  not 
decay  exponetially;  see  [6]  and  more  recently  [7].  In  such  cases  it  is  said  that  the  semigroup  does  not  satisfy  the 
spectrum-determined  growth  condition  [8].  The  determining  factor  in  the  examples  presented  in  [6]  and  [7]  can 
be  interpreted  as  the  accumulation  of  the  eigenvalues  of  A0  around  ±joo  in  the  form  of  Jordan  blocks  of  ever- 
increasing  size  (i.e.  as  the  eigenvalues  tend  to  ±joo  their  algebraic  multiplicity  increases  while  their  geometric 
multiplicity  stays  equal  to  one).  But  such  cases  are  ruled  out  when  one  deals  with  holomorphic  semigroups,  which 
is  the  reason  we  consider  these  semigroups  in  this  paper. 

Our  ultimate  aim  in  this  section  is  to  verify  exponential  stability.  By  Theorem  1,  in  order  to  prove  exponential 
stability  of  a  holomorphic  semigroup  with  infinitesimal  generator  A,  it  is  sufficient  to  show  that  S(4l)  C  C“. 
Hence,  in  the  first  part  of  this  section,  we  give  conditions  under  which  the  A  operators  described  by  (12)  generate 
holomorphic  semigroups.  In  the  second  part,  we  find  sufficient  conditions  which  guarantee  S(4l)  C  C“. 

Once  again,  the  setup  is  that  of  (12)  where  e  is  a  small  complex  scalar.  In  addition,  assume  that  Ao{k)  € 
is  diagonalizable  for  every  fc  G  M. 

Conditions  for  Sectorial  Infinitesimal  Generator 

To  find  conditions  under  which  A  in  (12)  will  define  a  holomorphic  semigroup  we  have  to  check  the  condition 
(18).  Since  the  Fourier  transform  preserves  norms,  (18)  is  equivalent  to  checking  \\{zl  —  A)~^\\  <  M/\z  —  a\  for 
z  belonging  to  some  sector  of  C.  This  involves  finding  the  inverse  of  the  operator  zl  —  A  and  then  calculating 
its  norm.  Such  a  computation  can  be  very  difficult  since  {zl  —  A)~^  has  the  form  (5)  [also  depicted  in  Figure  1 
(right)],  i.e.,  it  has  an  infinite  number  of  impulse  sheets.  On  the  other  hand,  finding  \\{zl  —  is  very  easy, 

due  to  the  diagonal  structure  of  A°.  Indeed  \\{zl  —  =  sup/^gj}  \\{zl  —  Ylo(A:))“^ ||. 

Thus  to  establish  conditions  for  A  to  be  sectorial,  we  again  use  perturbation  theory.  We  first  find  conditions 
under  which  A°  is  sectorial.  We  then  show  that  A  =  A°  eE  remains  sectorial  if  E  is  ‘weaker’  than  in  a 
certain  sense  we  will  describe,  and  if  e  is  small  enough. 

In  the  next  theorem  we  present  a  condition  for  a  spatially  invariant  A°  with  symbol  A^fk)  to  be  sectorial. 

Theorem  3:  Let  A{){k)  be  diagonalizable  for  every  fc  G  M,  and  let  U{k)  be  the  transformation  that  diagonalizes 
A(i{k),  i.e.,  4lo(fc)  =  U{k)K{k)  U~^{k)  with  A(fc)  diagonal.  Let  n{k)  :=  ||[/(fc)||  ||f7“^(fc)||  denote  the  condition 
number  of  U{k).  If  sup^gK  <  00,  and  for  every  fc  G  K  the  resolvent  set  p(Ao(fc))  contains  a  sector  of  the 
complex  plane  |  arg(z  —  a)  |  <  |  -f  7  with  7  >  0  and  a  G  M  both  independent  of  k,  then  A°  is  sectorial. 

Proof:  See  Appendix.  ■ 

This  theorem  has  a  particularly  simple  interpretation  when  Ao(-)  is  scalar.  In  this  case  K{U{k))  =  1  for  all 
A:  G  M.  Now  since  Ao(  • )  traces  a  curve  in  the  complex  plane,  by  Theorem  (3)  if  this  curve  stays  outside  some 
sector  I  arg(z  —  a)  |  <  |  -f  7,  7  >  0,  of  the  complex  plane  then  A°  is  sectorial. 
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The  following  theorem  is  from  [11].  It  uses  the  notion  of  relative  boundedness  [9]  of  one  unbounded  operator 
with  respect  to  another  unbounded  operator. 

Theorem  4:  Suppose  A°  S  ct,  M)  and  E  =  B°  F  C°  is  relatively  bounded  with  respect  to  A°  so  that 

\\E^\\  <  aUW  +  b\\A°i^\\,  (34) 

with  0  <  a  <  oo  and  0  <  6  |e|  <  1/(1  +  M).  Then  A  =  A°  +  eE  k  d,  sectorial  operator.  ■ 

This  theorem  says  that  if  A°  is  sectorial,  then  so  is  Yl  =  Yl°  +  e  if  is  weaker  than  A°  in  the  sense  of  (34) 
and  if  |e|  is  small  enough.  Notice  that  at  this  point,  condition  (34)  can  not  be  reduced  to  a  condition  in  terms 
of  Fourier  symbols  (i.e.  a  condition  that  can  be  checked  pointwise  in  k)  as  in  Theorem  3.  This  is  because  E  is 
not  spatially  invariant.  But  once  the  exact  form  of  the  operators  B°  and  C°  is  known,  (34)  can  be  simplified  to 
a  condition  on  the  symbols  of  A°,  B°  and  C°.  Let  us  clarify  this  statement  with  the  aid  of  an  example. 

Example  7:  Consider  the  periodic  PDF 

dt'4’  =  +  —  cip  +  e  dx  cos{flx)  Ip  +  u 

y  =  Ip- 

It  is  easy  to  see  that  A°  =  —(9^  +  —  c,  B°  =  dx  and  C°  =  6{x)  (the  identity  convolution  operator).  By 

formal  differentiation  we  have 

E  Ip  =  dx  cos{idx)  Ip  =  —  n  sin(nx)  t/i  +  cos{idx)  dx  tp ■ 

Using  the  triangle  inequality  and  ||  sin(na;)||  =  ||  cos(nx)||  =  1  we  have 

\\Ei;\\  <  \n\\m  +  WdxP^l  (35) 

Thus  we  have  effectively  ‘commuted  out’  the  (bounded)  spatially  periodic  operator  in  E,  and  are  left  with  only 
spatially  invariant  operators  on  the  right  of  (35).  Hence,  after  applying  a  Fourier  transfomation  to  the  right  side 
of  (35),  a  sufficient  condition  for  (34)  to  hold  is  that 

|H|  +  |fc|  <  a  +  &  +  c|,  k  G  M, 

which  can  be  shown  to  be  satisfied  for  large  enough  a  >  |H|  and  5  >  0.  ■ 

Remark  10:  The  above  example  makes  clear  the  notion  of  E  being  ‘weaker’  than  A°  that  we  mentioned  at 
the  beginning  of  this  subsection.  If  in  Example  7  we  had  B°  =  9^  and  C°  =  9/f  and  +  /i  =  5,  then  E  would 
contain  a  5*  order  derivative,  whereas  the  highest  order  of  dx  in  A°  is  4.  This  would  mean  that  (34)  can  not  be 
satisfied  for  any  choice  of  a  and  b.  m 

Conditions  for  Infinitesimal  Generator  to  have  Spectrum  in  C~ 

The  final  step  in  establishing  exponential  stability  is  to  show  that  S(A)  C  C“.  Unfortunately  it  is  in  general 
very  difficult  to  find  the  spectrum  of  an  infinite-dimensional  operator.  Thus  we  proceed  as  follows.  We  consider 
the  (block)  diagonal  operators  A°0,  9  G  [0,H).  This  allows  us  to  extend  Gersgorin-type  methods  (existing  for 
finite-dimensional  matrices)  to  find  bounds  on  the  location  of  S(^e),  A0  =  A°q  +  eSg.  In  turn,  we  use  this  to 
find  conditions  under  which  T,{A0)  C  C“,  and  thus  S(Yl)  C  C“. 

In  locating  the  spectrum  of  a  finite-dimensional  matrix  T  G  the  theory  of  Gersgorin  circles  [18]  provides 

us  with  a  region  of  the  complex  plane  that  contains  the  eigenvalues  of  T.  This  region  is  composed  of  the  union 
of  q  disks,  the  centers  of  which  are  the  diagonal  elements  of  T,  and  their  radii  depend  on  the  magnitude  of  the 
off-diagonal  elements  [18].  This  theory  has  also  been  extended  to  the  case  of  finite-dimensional  block  matrices 
(i.e.,  matrices  whose  elements  are  themselves  matrices)  in  [19].  Next,  we  further  extend  this  theory  to  include 
bi-infinite  (block)  matrices  Ag. 

For  the  operator  A  =  A°  +  e  E,  take  to  be  the  set  of  complex  numbers  z  that  satisfy 

CTmin(Z-f-  ^o(ft))  <  | e|  ( 1  (ft)  ||  -f  ||  (fc)  ||  )  , 

where  —  Ao(k))  denotes  the  smallest  singular  value  of  the  matrix  zl  —  Ao(k). 


(36) 
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Fig.  9.  Left:  The  regions  viewed  in  the  ‘complex-plane  X  spatial-frequency’  space  (the  disks  are  parallel  to  the  complex  plane). 

Center:  E(^9)  is  contained  inside  the  union  of  the  regions  Right:  The  bold  line  shows  E(j4“)  and  the  dotted  region  contains  E(A), 

A  =  A°  +  eE. 

Lemma  5:  For  every  9  we  have  T,{Ag)  C  &g,  where 

6,  =  U 

Proof:  See  Appendix.  ■ 

Example  8:  Let  us  consider  the  periodic  PDF  [5] 

dtf  =  +  —  cf  +  e  cos{flx)  doc  Ip  +  u 

y  =  (37) 

We  have 

Ao{k)  = -{k^  -  -  c,  B°{k)  =  l,  C°{k)=jk,  B{k)  =  l,  C{k)  =  l,  L  =  ^. 

From  (15),  Ai{k)  =  ^{k  -  fl),  A_i(fc)  =  +  fl),  and  thus  |jA_i(fc)||  +  ||Ai(A:)||  =  ^  (|fc  -  n|  +  |A:  +  n|). 

Hence  (36)  leads  to 

aUz/-Ao(fc))  =  \zI-Ao{k)\  <  M(|A:_H|  +  |fc  +  H|)  =  |  j^j  |  , 

which  means  that  the  set  &g  is  composed  of  the  union  of  disks  with  centers  at  Ao(0„)  and  radii  ^(|0n  —  H|  + 
\9n  +  H|).  This  is  nothing  but  an  extension  of  the  classical  Gersgorin  disks  to  bi-infinite  matrices.  Figure  9  (left 
&  center)  show  &g  in  the  complex-plane  x  spatial-frequency  space  and  in  C  respectively.^  ■ 

Remark  11:  The  set 

S,(M)  :=  {z  G  C  I  -  M)  <  s}  (38) 

=  {z  G  C  I  11(2/  —  M)ip\\  <  e  for  some  ||(/9||  =  1} 

=  {z  G  C  I  z  G  Sp(M  +  Z)  for  some  \\Z\\  <  e} 

is  called  the  e-pseudospectrum  of  the  matrix  M  [20].  Clearly  Yigi{M)  C  Se(M)  if  e'  <  e,  and  Se(M)  =  Sp(M) 
for  e  =  0.  The  pseudospectrum  is  composed  of  small  sets  around  the  eigenvalues  of  M.  For  instance  if  M  has 
simple  eigenvalues,  then  for  small  enough  values  of  e  the  pseudospectrum  consists  of  disjoint  compact  and  convex 
neighborhoods  of  each  eigenvalue  [21].  ■ 

^We  would  like  to  point  out  that  Figure  9  (left)  is  technically  incorrect;  once  the  spatially  invariant  system  is  perturbed  by  a  spatially 
periodic  perturbation  it  is  no  longer  spatially  invariant  and  thus  can  not  be  represented  by  a  Fourier  symbol.  Hence  its  spectrum  can  no  longer 
be  demonstrated  in  the  complex-plane  x  spatial-frequency  space.  Figure  9  (center)  demonstrates  the  correct  representation  of  the  Gersgorin 
disks  for  Ae- 
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Remark  12:  Comparing  (38)  and  the  definition  of  *8^  in  (36),  it  is  clear  that  Q5fc  =  Se(Ao(A:))  with  e  = 
|e|  (||y1_i(A:)||  +  ||7li(fc)||).  Thus  for  every  fc  S  M,  (36)  defines  a  closed  region  of  C  that  includes  the  eigenvalues 
of  A(i{k).  ■ 

We  now  employ  Lemma  5  to  determine  whether  S(y1)  resides  completely  inside  C“,  as  needed  to  assess  system 
stability.  Take  Dg  to  be  the  closed  disk  of  radius  e  and  center  at  the  origin,  and  18^  to  be  the  region  described 
by  (36).  Define  the  sum  of  sets  by  ifi  +  if2  =  {^  |  z  =  zi  +  Z2,  zi  G  iti,  2:1  G  ili}.  Also,  for  every  fc  G  M  let 
Amax(A)  represent  the  eigenvalue  of  Ao{k)  with  the  maximum  real  part,  and  let  K{k)  be  defined  as  in  Theorem  3. 
Theorem  6:  For  every  k,  is  contained  inside  Sp(Ao(fc))  +  S)r(fc)  with 

r{k)  =  \e\{\\A_,{k)\\  +  \\A,{k)\\)n{k). 

In  particular,  if  S(Gl°)  C  C“  and 

r{k)  <  |Re(An,ax(fc))|  + /?  (39) 

for  every  fc  G  M  and  some  P  <  0  independent  of  k,  then  S(A)  C  C“. 

Proof:  See  Appendix.  ■ 

Example  9:  Once  again  we  use  the  scalar  system  of  Example  8.  nfk)  =  1  since  A[){k)  is  scalar,  |Re(Aniax(fc))|  = 
|(fc^  —  +  c|,  and 

||A_i(fc)||  +  ||Ai(fc)||  =  l(|fc-0|  +  |fc  +  0|). 

Thus  condition  (39)  becomes 

\A  (^\k  -  n\  +  \k  +  n\)  <  \{k^->cY  +  c\  +  p. 

If  this  condition  is  satisfied  for  some  P  <  0,  the  dotted  region  in  Figure  9  (right)  will  remain  inside  C“  and  thus 

S(A)  c  C-.  ■ 

To  recap,  to  assess  exponential  stability  we  first  find  sufficient  conditions  on  A  such  that  it  belongs  to  the  class 
of  operators  for  which  the  spectrum-determined  growth  condition  holds.  These  are  conditions  under  which  A  is 
sectorial.  We  then  find  sufficient  conditions  for  A  to  have  C“  spectrum.  We  do  this  via  an  extension  of  Gersgorin 
circles  to  bi-infinite  (block)  matrices. 

VIII.  Conclusions  and  Future  Work 

We  use  perturbation  analysis  to  find  a  computationally  efficient  way  of  revealing  trends  in  the  Tf^-norm  of 
spatially  periodic  systems.  We  show  that  for  certain  classes  of  systems,  the  periodicity  can  be  chosen  so  as  to 
increase  the  Tf^-norm  or  induce  parametric  resonance.  An  application  of  this  would  be  in  mixing  problems.  It 
is  also  shown  that  the  Tf^-norm  can  be  made  to  decrease  for  an  appropriate  choice  of  the  frequency  of  the 
perturbation.  This  would  be  the  desired  scenario  in  the  design  of  the  body  of  aircraft.  We  demonstrate  that  for 
certain  scalar  systems,  the  value  of  the  spatial  period  that  achieves  the  desired  increase/decrease  of  the  Tf^-norm 
can  be  characterized  exactly  based  on  the  description  of  the  nominal  system. 

The  methods  presented  here  can  also  be  used  in  systems  with  many  spatial  directions.  For  example,  consider 
the  PDE 

ft  =  fyy  +  fxx  +  Clp  +  e  COs{flx)'tjj 

with  y  G  [—1,1]  and  a;  G  i?.  To  put  this  system  into  the  developed  framework  one  would  only  have  to  perform  a 
discrete  approximation  of  the  operator  dy  with  the  appropriate  boundary  conditions. 

We  also  study  the  problem  of  checking  the  exponential  stability  of  a  spatially  periodic  system.  We  do  this 
by  (i)  finding  conditions  under  which  the  A-operator  is  sectorial  (i.e.,  generates  a  holomorphic  semigroup)  and 
thus  satisfies  the  spectrum-determined  growth  condition,  and  (ii)  deriving  conditions  which  guarantee  that  A  has 
spectrum  contained  inside  the  open  left-half  of  the  complex  plane. 

Future  research  in  this  direction  would  include  an  exact  (analytical)  characterization  of  the  frequencies  for 
which  the  Tf^-norm  is  most  increased/decreased  for  the  general  case  of  matrix-valued  A(-).  The  perturbation 
methods  presented  here  could  also  be  generalized  to  bi-infinite  Sylvester  equations  which  arise  frequently  in  fluids 
problems. 
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IX.  Appendix 

Proof  of  Theorem  3 

It  is  shown  in  [22]  that  a  sufficient  condition  for  A°  to  be  sectorial  is  that  p{A°)  contain  some  right  half  plane 
{z  G  C  I  Re(z)  >  p},  and 

\\z{zl  -  A°)-^\\  <  M  for  Re(z)  >  p, 


for  some  /x  >  0  and  M  >  1. 

Now  since  Ao{k)  G  has  simple  eigenvalues  for  every  k,  there  exists  a  matrix  U{k)  such  that  A[){k)  = 

U{k)  A(fc)  U~^{k)  with  A(fc)  a  diagonal  matrix.  Let  Xi{k),  i  =  1,  •  •  •  ,q  denote  the  diagonal  elements  of  A(fc). 
Clearly  these  are  also  the  eigenvalues  of  Ao(A:).  Thus  we  have 

\\z{zI-A°)-^  <  snp(\\z{zI-Ao{k))-^) 

kGR  ^  ^ 

<  sup(||C/(fc)||  \\U-Hk)\\  \\z{zl  -  A{k))-^) 

kGR  ^  '' 

=  dist[z,llpU(*))]) 

/  \z\  \ 

f:  ^max  sup  (  .  „  .  .  TTTvT  j  ’ 

km  Vdist[z,Sp(Ao(A:))]/ 

where  Kmax  :=  supj.gjj  K{k). 

Let  us  now  choose  the  positive  scalar  M'  =  (1  +  Kmax)Af,  M  >  1,  and  consider  for  a  given  k  the  region  of 
the  complex  plane  where 

.  _ N _  > 

dist[z,Sp(Ao(Ai))]  -  ■ 


This  region  (which  contains  the  eigenvalues  Xi{k))  is  contained  inside  the  union  of  the  circles 

i  =!,■■■  ,q, 


i  =  l,  •••,<?•  (Al) 

Notice  that  (Al)  describes  circles  whose  radii  increase  like  \Xi{k)\/M,  M  >  1,  as  their  centers  Xi{k)  become 
distant  from  the  origin.  Clearly  a  sufficient  condition  for  these  circles  to  belong  to  some  open  half  plane  {z  G 
C  I  Re(z)  <  /x}  for  all  fc  G  M  and  large  enough  M  is  that  Sp(Ao(A:)),  fc  G  M,  reside  outside  some  sector 
I  arg(z  —  a)  I  <  §  +  7,  7  >  0,  of  the  complex  plane. 

Finally,  if  the  circles  (Al)  are  contained  in  some  open  half  plane  {z  G  C  |  Re(z)  <  p}  for  all  fc  G  M,  then  for 
Re(z)  >  ^,  z  G  p(Ao(A:))  and  we  have 


t^max  I  ,  /,N|  ^  Af  , 

|z  -  Xi(k)\ 

which  are  themselves  contained  inside  the  larger  circles 

|A*(fc)| 


\z-X^{k)\  < 


M 


^max 


sup 

kGR 


( _ y ) 

Vdist[z,  Sp(Ao(fc))]/ 


and  thus  ||z(z/  —  A°)  ^\\  <  M  for  Re(z)  >  /x. 


<  M 


Proof  of  Lemma  5 

We  use  n^vTIIiv  to  denote  the  (2N  +  1)  x  (2N  +  1)  [block]  truncation  of  an  operator  T  on  (f,  where  IIjv 
is  the  projection  defined  by 

center 

diag|--  -  ••• 


2N+1  times 
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where  I  is  the  identity  matrix.  Notice  that  Hp^TUn  is  still  an  operator  on  it  made  from  the  bi-infinite 
T  by  replacing  all  entries  outside  the  center  (2A^-|- 1)  x  (2iV-|-l)  block  with  zeros.  We  now  form  the  finite¬ 
dimensional  matrix  n7v.4enArL  „  by  restricting  flTv^ellAr  to  the  finite-dimensional  space  If^vf^.  Clearly 
Hn  Ae  n^v  ri  ei  has  pure  point  spectrum.  Hence  using  a  generalized  form  of  the  Gersgorin  Circle  Theorem  [19] 

1  i  c 

for  finite-dimensional  (block)  matrices,  we  conclude  that 

S(nAr c  [J  Q3e„  C  [J 

|n|<Al  neZ 

where  Q3e„  are  regions  of  C  defined  by  (36).  Since  this  holds  for  all  TV  >  0,  we  have  S(-4e)  C  ©g. 

Proof  of  Theorem  6 

If  U{k)  diagonalizes  Ao{k),  Ao{k)  =  U{k)A{k)  U~^{k),  and  K{k)  =  ||{7(fc)||  ||C/“^(fc)||  denotes  the  condition 
number  of  U{k),  then  from  [23]  the  pseudospectrum  of  Ao{k)  satisfies 

^p(^o(^))  +  ®e  ^  '^s{Ao{k))  C  Sp(ylo(fc))  -f  (A2) 

for  all  £  >  0.  Thus  the  first  statement  of  the  Theorem  follows  immediately  from  (A2)  and  the  fact  that  = 
Se(Ao(fc))  with  £  =  |e|  (||£l_i(fc)||  +  ||Ai(fc)||)  [see  Remark  12].  To  prove  the  second  statement,  let  denote 
all  complex  numbers  with  real  part  less  than  /3  G  M.  It  follows  from  S(£l°)  C  C“  that  C  C“  for  every 

9.  If  (39)  holds  then 

C  Sp(£lo(0ri))  +  ®r(e„)  C 

for  every  n  G  Z,  and  from  Lemma  5  we  have  ^^{Ae)  C  C  C^,  for  some  /?</?'<  0  and  every 

6».  Thus  S(A)  C  C". 


Interpretation  of  Imaginary  States 

It  was  shown  above  that  one  can  decrease  the  Tf^-norm  of  certain  systems  by  choosing  the  perturbation 
amplitude  Ai  to  be  purely  imaginary  (or,  in  general,  skew-symmetric).  This  would  yield  a  perturbed  system  that 
can  in  general  have  states  with  nonzero  imaginary  parts.  One  could  then  ask  the  physical  interpretation  of  such 
a  system. 

For  any  operator  Ap  and  function  ip  one  can  write 

Ap  =  Ar  +  jAi,  Ar,  Ai  G 


Then  the  system  equations  can  be  written  as 

dtifr  +  jfi)  =  {Ar  +  jAi){lpr  +  jfi) 

t 

—  Arlpr  Ailpi 
—  Arlpi  Ailpr 

t 

r  1  r 

dt 

Clearly  the  state  dimension  is  twice  that  of  the  original  system  with  imaginary  coefficients,  but  now  [^''j  G 
Let  us  give  a  simple  example.  Assume  the  heat  equation,  with  Ai  =  je,  e  G  M,  i.e.  Ap  =  —  c  +  jecos{^lx). 

Then  Ar  =  Aq  =  —  c,  Ai  =  Ai  =  jecos{^lx),  and  thus 


Ipr 

Ar  —Ai 

Ipr 

'lpi_ 

Ai  Ar 

fi  _ 

dt 


'4^r 

dl-c 

— ecos(Ha;) 

'4^r 

e  cos  (Ox) 

dl-c 

dt-c  0 


-I-  £cos(Ha;) 


0  -1 
1  0 


0  dl  —  c 

which  describes  two  identical  systems  coupled  through  the  periodic  perturbation. 
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